function [Kr2_1,Kr2_2,Kr2_3,Kr2_4] = lpcvd_Kr2_sincos(L,mode)
Kr2_1 = ones(mode+1,mode+1)/L^2;
% Kr2_1(1,:) = 0;
% Kr2_1(:,1) = 0;

Kr2_2 = Kr2_1;
Kr2_2(1,1) = 0.0;

Kr2_3 = Kr2_1;
% Kr2_3(1,:) = 0;

Kr2_4 = Kr2_1;
% Kr2_4(:,1) = 0;
% for m = 0:mode
%     for n = 0:mode
%         if m==0 || n==0
%             Kr2_1(m+1,n+1) = 0;
%         else
%             Kr2_1(m+1,n+1) = 1/4;
%         end
%         
%         if m==0 && n==0
%             Kr2_2(m+1,n+1) = 0;
%         elseif m==0 || n==0
%             Kr2_2(m+1,n+1) = 1/2;
%         elseif m~=0 && n~=0
%             Kr2_2(m+1,n+1) = 1/4;
%         end
%         
%         if m == 0
%             Kr2_3(m+1,n+1) = 0;
%         elseif m~=0 && n==0
%             Kr2_3(m+1,n+1) = 1/2;
%         elseif m~=0 && n~=0
%             Kr2_3(m+1,n+1) = 1/4;
%         end
%         
%         if n == 0
%             Kr2_4(m+1,n+1) = 0;
%         elseif n~=0 && m==0
%             Kr2_4(m+1,n+1) = 1/2;
%         elseif n~=0 && m~=0
%             Kr2_4(m+1,n+1) = 1/4;
%         end
%     end
% end